Anomalous lattice response at the Mott Transition in a Quasi-2D Organic Conductor 
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Discontinuous changes of the lattice parameters at the Mott metal-insulator transition are detected 
by high-resolution dilatometry on deuterated crystals of the layered organic conductor re-(BEDT- 
TTF)2Cu[N(CN)2]Br. The uniaxial expansivities uncover a striking and unexpected anisotropy, 
notably a zero-effect along the in-plane c-axis along which the electronic interactions are relatively 
strong. A huge thermal expansion anomaly is observed near the end-point of the first-order transition 
line enabling to explore the critical behavior with very high sensitivity. The analysis yields critical 
fluctuations with an exponent a ~ 0.8 ± 0.15 at odds with the novel criticality recently proposed 
for these materials [Kagawa et al, Nature 436, 534 (2005)]. Our data suggest an intricate role of 
the lattice degrees of freedom in the Mott transition for the present materials. 
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The Mott metal-insulator (MI) transition has been the 
subject of intensive research for many years, see e.g. [lj 
for a review. Materials intensively discussed in this con- 
text include transition metal oxides, notably Cr-doped 
V2O3, and, recently organic k-(BEDT-TTF) 2 X charge- 
transfer salts @, H |j, Here BEDT-TTF (or simply 
ET) denotes bis(ethylenedithio)tetrathiafulvalene and X 
a monovalent anion. For the latter substances, pressure 
studies revealed a first-order metal-insulator transition 
line T MI (P) @, 0, 0, & @, indicative of a bandwidth- 
controlled Mott transition 0, and suggest a second- 
order critical endpoint (P ,T ) 0, 0, 0, [f| with remark- 
able properties. Particularly striking was the observation 
of a pronounced softening of the C22 elastic mode UJt. Al- 
though acoustic and lattice anomalies are expected|9l. Iioj 
at the Mott transition in response to the softening of the 
electronic degrees of freedom, the actual role of the lattice 
for the Mott transition in real materials remains illusive. 
In addition, an unconventional Mott criticality was pro- 
posed for the present organic salts [ll[ and attributed to 
their quasi-twodimcnsional (quasi-2D) electronic charac- 
ter. 

In this Letter we report, for the first time, the direct 
observation of lattice anomalies at the Mott transition in 
a k-(ET)2X organic conductor and explore, via a sensi- 
tive thermodynamic probe, the criticality near (P ,T ). 

For the thermal expansion measurements, a high- 
resolution capacitive dilatometer (built after was 
used, enabling the detection of length changes Al > 
1CP 2 A. Owing to the experimental difficulties posed 
by accomplishing high-resolution dilatometric measure- 
ments under variable pressure, use was made of the pos- 
sibility of applying chemical pressure. To this end, single 
crystals of K-(d8-ET) 2 Cu[N(CN) 2 ]Br were synthesized 
with deuterium atoms replacing the protons in the ter- 
minal ethylene groups. These fully deuterated salts, re- 
ferred to as d8-Br in the following, are known to be sit- 
uated very close to the MI transition First, deuter- 



ated (98%) ET molecules were prepared according to 



14, 15] using multiple rccrystallization for the intermedi- 



ate steps. Next, single crystals were synthesized along an 
alternative preparation route described recently for the 
protonated variant h8-Br [l|| . The grade of deuteration 
was checked by infrared reflection spectroscopy both on 
the deuterated ET material [13] as well as on the d8- 



Br single crystals [18j, and found to be at least 98%. For 
the present study, crystals of two independently prepared 
batches were used: crystal #1 (batch A2907) and #3 
(A2995). The crystals have the shape of flat distorted 
hexagons with dimensions of about 1 x 1 x 0.4 mm 3 . 
The pressure exerted on the crystal by the dilatometer 
ranges from 1 to 6 bar. A preliminary account on a sec- 
ond crystal from batch A2907 was given in (19^ . The 
resistivity was studied by employing a standard four- 
terminal ac-technique. All measurements, unless stated 
otherwise, were carried out after cooling through the 
glass transition at T g ~ 77 K with a very low rate of 
-3K/h (thermal expansion) and -6K/h (resistivity) to 
rule out cooling-rate dependent effects, see [2oj ] . 

The interlayer resistivity p± for crystal #1 is shown in 
the lower part of Fig.[T] Upon cooling, p± passes over a 
maximum around 45 K, then rapidly drops and flattens 
around 30 K. The resistivity remains metallic down to 
about 20 K, below which the slope sharply increases (cf. 
upper inset in Fig.[l| indicating the transition into an 
insulating state. A similar 01 was found for crystal #3 
and the crystal studied in [19] including the vanishing 
of p± below about 11.5 K. A zero resistivity accompa- 
nied by a tiny signature in the oti data is consistent with 
percolative superconductivity in a minor metallic phase 
coexisting with an antiferromagnetic/insulating ground 
state for d8-Br [21(, cf. the phase diagram in Fig. [31 

The main features in the resistivity have their clear 
correspondence in the coefficient of thermal expansion, 
a = l^dl/dT, also shown in Fig. Q] along the a-axis of 
crystal #1. The flattening of p± is accompanied by a 
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FIG. 1: (Color online) In-plane a-axis expansivity, a a (full cir- 
cles), (left scale) and interlayer resistivity, p±, (right scale) for 
K-(d8-ET) 2 Cu[N(CN) 2 ]Br crystal #1 together with a a data 
for crystal #3 (open circles). Upper inset: blow-up of low-T 
p± (T) data on the same T scale as used in the main panel. 



huge peak in a(T) centered at a temperature referred 
to as T p in the following around 30 K. As will be dis- 
cussed below, this effect can be assigned to a second-order 
phase transition. Upon further cooling, a a (T) reveals an 
even bigger negative peak indicating yet another phase 
transition. The accompanying change in pj_ from metal- 
lic to insulating behavior suggests this peak to be due 
to the MI transition. This is consistent with measure- 
ments under magnetic fields up to 10 T (not shown) leav- 
ing the peak position unaffected. A very similar a a (T) 
behavior is observed for #3, although with slightly re- 
duced (<~ 20%) peak anomalies and minor shifts in T p 
and T MI , cf. Fig.Q] More insight into the character of 
the transitions can be gained by looking at the relative 
length changes Ak(T)/k = (k(T) - Z. ( (300 K))/k(300 K), 
(i = a, b, c) shown in Fig.[2]for crystal #1. The dominant 
effects occur along the in-plane a-axis, i.e. parallel to the 
anion chains. Here a pronounced s-shaped anomaly is 
revealed at T p which lacks any sign of hysteresis upon 
cooling and warming - generic features of a second-order 
phase transition with strong fluctuations. On further 
cooling through Tmi, the a-axis shows a rapid increase 
of about Aa/a = 3.5 • 10~ 4 within a narrow temperature 
range, indicative of a slightly broadened first-order tran- 
sition. The observation of a small but significant hystere- 




T(K) 

FIG. 2: Relative length changes for fi>(d8-ET) 2 Cu[N(CN) 2 ]Br 
crystal #1 along the in-plane a- and c-, and out-of-plane 
b-axis. The data have been offset for clarity. Lower inset 
shows hysteretic behavior in Aa/a at Tmi measured at very 
low sweeping rates of ±1.5 K/h. Upper inset depicts the 2D 
triangular-lattice dimer model with transfer integrals t and t' . 

sis of about 0.4 K (cf. lower inset Fig. [2]), which complies 
with the hysteresis in p±(T) (upper inset Fig.[T]), confirms 
the first-order character. The corresponding anomalies 
along the 6-axis are less strongly pronounced. Surpris- 
ingly, for the second in-plane c-axis, anomalous behavior 
in Al/l can be discerned neither at T p nor at Tmi- The 
same anisotropy was found for the second crystal of batch 
A2907 studied in [l9| (not shown), on which all three uni- 
axial expansion coefficients had been determined. 

Figure [2] reveals that the anomalies at T p and Tmi are 
correlated in size, albeit with reversed sign, suggesting 
that they arc intimately related to each other. In ad- 
dition, the data disclose a striking in-plane anisotropy. 
Given the quasi-2D electronic structure as shown in the 
upper inset of Fig. [21 characterized by dimers on an 
anisotropic triangular lattice the latter is a very re- 
markable and unexpected result: The dominant response 
in the a-axis, along which no direct dimer-dimer overlap 
exists, means that the diagonal electronic interactions 
along the c ± a directions, t, have to be involved in this 
process. Since these interactions have a strong compo- 
nent also along the c-axis, which is likely to be even softer 
than the anion-chain a-axis [22], a significant c-axis re- 
sponse would be expected at Tmi- A zero-effect along 
the c-axis is even more amazing as there is a relatively 
strong direct dimer-dimer interaction t' along this axis, 
cf. upper inset of Fig. [2] Thus, to account for a zero c-axis 
response within a 2D electronic model would imply an ac- 
cidental cancellation of counteracting effects associated 
with t and t' , which seems very unlikely. Furthermore, it 
is not obvious how these in-plane interactions may cause 
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the comparatively strong effect in the interlayer 6-axis, 
along which the lattice is expected to be even more stiff 
|22l ]. These observations suggest that a coupling of the 
7r-electrons to other degrees of freedom has to be taken 
into account to understand the MI transition here. 

Before discussing further implications of our observa- 
tions, the Mi-transition temperature Tmi is determined. 
As Fig. [2] demonstrates, the transition is not very sharp 
but rather spans a range of several Kelvin - an effect 
which is very similar for both crystals studied here and 
the one explored in A broadening of signatures 

in T-dependent measurements, as opposed to isothermal 
pressure sweeps, would be naturally expected given the 
steepness of Tmi(P), cf. Fig.[3l However, the width of 
about 5.6 K of the Al a /l a jump (10-90%), which trans- 
forms into a pressure interval of about 2 MPa employing 
a slope dT MI /dP = -2.7 ± 0.1 K/MPa around 14K (cf. 
Fig. [3]), is even smaller than the transition range seen in 
acoustic measurements as a function of pressure [H (cf. 
hatched area in Fig. [3]), but is comparable with a width 
of about 1.4 MPa as read off the resistivity data in 
These smearing effects have been attributed to a region 
of coexistence between insulating and metallic phases Q , 
as indeed observed via real-space imaging [24|. For lack 
of a well-founded procedure to treat the broadened tran- 
sitions, the position of the a(T) minimum is chosen as 
the thermodynamic transition temperature. Employing 
literature results on Tmi(P) (Hj, the so-derived values 
of T M i = (13.5 ± 0.8) K, (#1) and (14.1 ± 0.8) K, (#3) 
can be used to pinpoint the position of the present d8-Br 
crystals on the pressure axis in Fig. [3] Within the uncer- 
tainties implied in this procedure, the crystals are located 
very close to the critical pressure Pn as determined by the 
various pressure studies 0, 0, 0, Hfl • The significance of 
this finding is twofold. First, it demonstrates that the 
anomaly at T p reflects the lattice response at (Po,To). 
Second, as this point is part of the Tmi line, it provides 
a natural explanation for the intimate interrelation of the 
anomalies at T p and Tmi inferred from Fig. [2] 

The huge anomaly at T p ~ To, exceeding the back- 
ground by a factor 3-4, enables us to explore the critical- 
ity at (Po,To) with extraordinarily high sensitivity. To 
this end, the phase transition anomaly in a a (T), shown 
for the crystals #1 and #3 in Fig. [4] on expanded scales, 
is analyzed in terms of a power-law behavior in the vari- 
able t = (T — To) /To. This approach is based on the 
proportionality of a(T) to the specific heat, C(T), im- 
plying that the same scalin g la ws apply at T , as verified 
by various groups, see, e.g. 2f| 27j |. 

The data sets in Fig.Q] reveal a steep increase in the 
slope of a(T) on the outer flanks of the maximum. Closer 
to the center of the peak, however, the slope is reduced 
giving rise to a rounded maximum. Such broadening ef- 
fects are generally encountered in the immediate vicin- 
ity of the transition and attributed to sample inhomo- 
geneities. The rounding over a considerable temperature 
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FIG. 3: (Color online) Phase diagram of k-(ET)2X including 
the Neel transition at Tjv (dashed line) and the superconduct- 
ing (SC) transition at T c (dotted line) from [2|. Thin solid 
lines, delimiting the hatched area, denote positions of acoustic 
anomalies Q. The middle position (thick solid line) is used 
here as Tmi{P) [H|]. Closed symbols refer to anomalies at 
Tmi and T p ~ To of the d8-Br crystals #1 and #3 (same po- 
sitions on the scale of the figure), while open symbols denote 
literature results for (P ,T ): (o) 0], (v) 1, (O) Q, (A) 
PJ. Vertical line indicates the T-sweeps performed here. 
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FIG. 4: (Color online) Expansivity along the a-axis for d8-Br 
crystals #1 (o) and #3 (v) near To. Solid (dashed) lines are 
fits as described in the text for a — 0.8 with (without for #1) 
a Gaussian distribution of To . Straight lines show background 
contributions implied in fits for a values given in the figure. 
Star marks universal background point discussed in the text. 



range here demands particular attention. For the de- 
scription of the data in the range 24 - 36 K, the function 
a = 4^|i| _a + B + E ■ T was used. This function con- 
tains the singular contribution with the amplitudes A + 
and A~ for t > and t < 0, respectively, and a linear 
term. The latter comes primarily from the phonons but 
can also include a small non-singular electronic contri- 
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bution. The smearing of the transition is accounted for 
by a Gaussian distribution for To, G(Tq,Tq, 5Tq), cen- 
tered at T with a width ±5T . Applying the function 
J [4f | (T - T )/T 1 - & + B + E ■ T] ■ G(f , T ,ST ) ■ dT , the 
data sets of crystal #1 and #3 were fitted simultaneously 
using the same exponent a, the same ratio A + /A~, and 
an identical background for both crystals. A constraint 
for the background contribution can be derived by com- 
paring the data in Fig.|H with those of the d8-Br in 
and h8-Br in ■ All data sets intersect at a single point 
T ~ 21 K, a ~ 16 • 10 _6 K _1 , irrespective of the presence 
and size of the critical contribution at To, indicating that 
this point reflects the pure background. Thus a mean- 
ingful background should extrapolate to this universal 
point. A good fit to both data sets, also satisfying this 
background constraint, is obtained for a — 0.8, A + /A~ 
= 0.79, and T = 30.1 K, 6T = 1.59 K for #1 and T = 
29.6 K, 5T = 1.74 K for #3, cf. Fig.g] We stress that 
a values in the range 0.65 - 0.95, with small changes in 
the other parameters accordingly, result in fits of similar 
quality and still comply with the background constraint. 
In contrast, the residual of the fit increases substantially 
upon decreasing a to well below 0.65. This is accompa- 
nied by a suppression of the background to even negative 
values for a < 0.5, clearly incompatible with the back- 
ground constraint, see, e.g., the background implied in 
the fits for a = 0.6 and 0.5 in Fig.[U As clearly indicated 
by these simultaneous fits, and confirmed by independent 
fits to the individual data sets for crystals #1 and #3, 
a large positive a value is the only possible, physically 
meaningful description of the expansivity data. 

The critical exponent derived here of a ~ 0.8 ± 0.15 
is much larger than those of known universality classes 
with -0.12 < a < 0.14 and the mean-field value a — 
observed at the Mott critical endpoint of Cr-doped V2O3 
[lH . In particular, it greatly conflicts with the criticality 
reported in [HI for pressurized X = Cu[N(CN) 2 ]Cl. Em- 
ploying the exponent identity a+2(3+j = 2 [29j the expo- 
nents found there of (S, 0, 7) rs (2, 1, 1) give a = -1. The 
reason for this discrepancy is unclear but might be re- 
lated to the significant broadening effects [IB] , which have 

The exponent 
= 0.5 expected 
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not been included in the analysis in 
found here, however, is rather close to a 
for a tricritical point 31| . Such a scenario would imply a 
symmetry breaking associated with Tmi for which no evi- 
dence has yet been supplied. Interestingly, an even larger 



exponent a = 0.93 was reported for Lao.7Cao.3Mn03 [271 ] . 
also characterized by a strong electron-phonon coupling, 
showing a similar a(T) anomaly as the one observed here. 

In summary, high-resolution dilatometry on deuter- 
ated K-(ET) 2 Cu[N(CN) 2 ]Br crystals reveals discontinu- 
ous changes of the lattice parameters at the Mott transi- 
tion. The data disclose a striking anisotropy unlikely to 
be captured by a 2D purely electronic model. An analysis 
of the huge thermal expansion anomaly at the end-point 



of the first-order Tmi hue yields a critical exponent a ~ 
0.8 ± 0.15, markedly different from the criticality derived 
from transport measurements [H| . The unusually large a 
value together with the anomalous anisotropy of the lat- 
tice effects at Tmi suggest an intricate role of the lattice 
in the Mott transition for the present materials. 
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Brazilian Research Foundation CAPES and the DAAD. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 
[12] 

[13] 
[14] 
[15] 

[16] 

[17] 
[18] 

[19] 



[20] 
[21] 
[22] 



[23] 
[24] 
[25] 



M. 

S. 

P. 

D. 

F. 

H. 

K. 

H. 



F. 
R. 
A. 
K. 
M. 



[26] 
[27] 
[28] 
[29] 
[30] 



[31] 



Imada et al., Rev. Mod. Phys. 70, 1039 (1998). 
Lefebvre et al, Phys. Rev. Lett. 85, 5420 (2000). 
Limelette et al, Phys. Rev. Lett. 91, 016401 (2003). 
Fournier et al, Phys. Rev. Lett. 90, 127002 (2003). 
Kagawa, et al, Phys. Rev. B 69, 064511 (2004). 
Ito et al, J. Phys. Soc. Jpn. 9, 2987 (1996). 
Kanoda, Hyperfine Interact. 104, 235 (1997). 
Kino et al, J. Phys. Soc. Jpn. 64, 2726 (1995). 
S.R. Hassan et al, Phys. Rev. Lett. 94, 036402 (2005). 
J. Merino et al, Phys. Rev. B 62, 16442 (2000). 
Kagawa et al, Nature 436, 534 (2005). 
Pott et al, J. Phys. E 16, 444 (1983). 
Kawamoto et al., Phys. Rev. B 55, 14140 (1997). 
Hartke et al, Chem. Ber. 113, 1898 (1980). 
Mizuno et al, J. Chem. Sci. Commun. 1978, 18 
(1978). 

Ch. Strack et al, Phys. Rev. B 72, 054511 (2005). 

S. Gartner et al., Synth. Met. 44, 227 (1991). 

E. Griefihaber, Dissertation, University Stuttgart (2000), 

unpublished. 

M. Lang et al, Proceedings of the 8th Intern. Conf. on 
Materials and Mechanisms of Superconductivity - High 
Temperature Superconductors, Physica C, in press. 
J. Miiller et al, Phys. Rev. B 65,' 144521 (2002). 
K. Miyagawa et al, Phys. Rev. Lett. 89, 017003 (2002). 
For the related X = Cu(NCS)2 salt, the uniaxial com- 
pressibilities ki are strongly anisotropic with ki : k-z : 
kz = 1 : 0.53 : 0.17, and i — 1 and 2 the in-plane axis per- 
pendicular and parallel to the anion chains, respectively, 
and i = 3 along the long axis of the ET molecules. 
D. Chasseau et al, Synth. Met. 41 - 43, 2039 (1991). 
T. Sasaki et al, Phys. Rev. Lett. 92, 227001 (2005). 
Due to the finite range of phase coexistence around Tmi, 
the comparison of the thermodynamically determined 
Tmi values here with those derived from transport ex- 
periment @, 0] may cause some ambiguity. We therefore 
refer to Tmi(P) as the mean value (thick solid line in 
Fig. [3} of acoustic anomalies Q. 
V. Pasler et al., Phys. Rev. Lett. 81, 1094 (1998). 
J. A. Souza et al, Phys. Rev. Lett. 94, 207209 (2005). 
P. Limelette et al, Science 302, 89 (2003). 
L.P. Kadanoff et al., Rev. Mod. Phys. 39, 395 (1967). 
In view of the huge dTo/dP of approximately 2-3K/MPa, 
cf. Fig. |3] weak strain fields due to minor sample inhomo- 
geneities can give rise to a broad distribution of To values 
as observed here. 

K. Huang, Statistical Mechanics (Wiley, New York, 1987) 
2 nd ed., Chap. 17.6. 



